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Abstract 

In this paper we analyze different contributions to the magnetoresistance of 
magnetic tunneling junctions at low voltages. A substantial fraction of the 
resistance drop with voltage can be ascribed to variations of the density of 
states and the barrier transmission with the bias. However, we found that 
the anomaly observed at zero bias and the magnetoresistance behavior at very 
small voltages, point to the contribution of inelastic magnon-assisted tunnel- 
ing. The latter is described by a transfer parameter T J , which is one or two 
orders of magnitude smaller than T d , the direct transmission for elastic cur- 
rents. Our theory is in excellent agreement with experimental data, yielding 
estimated values of T J which are of the order of T d /T J ~ 40. 
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I. INTRODUCTION 



Recently, the interest on the phenomena of Giant Magnetoresistance (GMR) in magnetic 
tunnel junctions (MTJ) has grown significantly due to potential applications in magnetore- 
sistive reading heads, magnetic field sensors, nonvolatile magnetic random access memories, 
and many others 1-5 . The effect is based on the spin dependent scattering mechanisms 
proposed in the early papers by Cabrera and Falicov 6 , which lead in MTJ's, to a strong 
dependence of the conductance on the magnetic polarization 7 . Typically, the GMR effect 
found in MTJ's is of the order of 25 — 30% 8 ' 9 , and points to a large ratio of the densities of 
states for majority (M) and minority (m) electrons at the Fermi level (E F ) 

N M (E F ) 

/ « 2.0 - 2.5 . 

N m (E F ) 

As usual in MR experiments, one compares the resistances for the cases where the magne- 
tizations at the electrodes are anti-parallel (AP) and parallel (P). In several experiments 
reported in the literature (see for example Ref. 1, 2, 8, 9), the junction resistance drops sig- 
nificantly with the applied voltages, with a sharp peak at zero bias (zero-bias anomaly). This 
bias dependence shows a rapid initial decrease up to voltages of the order of V ~ 100 mV, 
then slows down but continues decreasing with voltages, up to 60% of the peak value at 
500 mV in some cases 9 . Many attempts to explain the above behavior have been done over 
the last years 1 ' 2 ' 8 ' 10 , but a complete theory which includes all the observed features is still 
lacking. 

In Ref. 8, scattering from magnons at the electrode-insulator interface has been proposed 
as the mechanism for randomizing the tunneling process and opening the spin-flip channels 
that reduce the MR. While this process may explain the MR behavior in the vicinity of 
zero-bias (voltages smaller than 40 — 100 mV), estimations of magnon scattering cross 
sections show that the effect is too small to account for the sharp drop in resistance observed 
in the whole range of 500 mV. In fact, inelastic-electron tunneling spectroscopy (IETS) 
measurements at low temperature showed peaks which can unambiguously be associated 
with one-magnon spectra at very small voltages (from 12 to 20mV, with tails up to 40mV, 
and maximum magnon energy not larger than 100 meV) 1 . To go beyond this limit will imply 
multimagnon processes, which are negligible at low temperature. This way, the electron- 
magnon coupling constant coming from Ref. 8 is by sure considerable overestimated. 

The above explanation 8 has been challenged in Ref. 10, where it is shown that the 
experimental data can be understood in terms of elastic tunneling currents which conserve 
spin, by considering effects not taken into account in Ref. 8. Those include the lowering of 
the effective barrier height with the applied voltage, as in the classical Simmons' theory 11 , 
and most important, variations of the densities of states with the bias at both magnetic 
electrodes. The latter is a relevant question, since experiments probe depths of the order of 
0.5 eV from the Fermi surface. The simple calculation developed in Ref. 10 models the band 
structure with free electron-like densities of states, since the tunneling current is dominated 
by the s-electron contribution. This approach yields a zero bias anomaly which depends on 
the band structure, and a variation of the MR which has the right order of magnitude for 
the whole range of 500 meV. The above discussion and other experimental results primarily 
exhibit that the density of states dependence on the applied voltage plays an important 
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role 12 ' 13 . However, fine details of experiments at very small voltages are difficult to fit. One 
may adopt here a pragmatic procedure, with a more intricate band structure and more free 
parameters to improve the fitting 2 . 

In this paper we take a different stand, motivated by results from IETS experiments 1 , 
which show that inelastic scattering do participate in the phenomenon at very small voltages. 
Also, MR experiments 8 ' 9 show clearly a different behavior with applied voltage in the same 
small bias region (up to 100 mV). A complete theory then should include: 

i) magnon assisted tunneling effects, with maximum magnon energies of the order of 
~ 100 meV. At low temperature, electrons from the electrodes, accelerated by the 
applied voltage, excite magnons at the interface. At low temperature, only magnon- 
emission processes should be considered; 

ii) variation with voltages of the densities of states for the different spin bands in the 
ferromagnets. Here, we will follow closely the approach of Ref. 10, with a simple 
picture of the band structure. This is motivated by the discussions given in Ref. 14 
and 10 over the polarization of the tunneling current. We assume here that the latter 
is mainly of s-character; 

Hi) lowering of the effective barrier height with the applied voltage. This effects, as shown 
in Ref. 11, yields to a parabolic dependence of the resistance with the bias. It does 
not contribute to the zero-bias anomaly, but it is always present and should dominate 
the behavior at large voltages. 

The above program will be developed in the present contribution. The content of this 
paper can be described as follows: In the next Section, we formulate the theoretical basis 
for analyzing tunneling currents, discussing the transfer Hamiltonian which includes all 
the above mentioned ingredients. In Section 3, we solely analyze density of states effects, 
considering elastic tunneling processes. Some analytical expressions are shown. In Section 4, 
we include contributions from inelastic magnon- assisted processes to the tunneling current, 
and finally, in the last Section, a few conclusions and remarks are added. 

II. THEORETICAL PRELIMINARIES 

To give a description of the MR and the resistance in the MTJ, we will use the transfer 
Hamiltonian method 15 . The junction is composed by two ferromagnetic electrodes separated 
by a thin oxide film which represents a potential barrier due to the fact that the Fermi 
levels of the ferromagnetic layers are situated in the gap region of the oxide film. We have 
considered the s-band electrons as free particles (plane- waves) , being responsible for the 
dominant contribution to the tunneling process. The d electrons, which are more localized, 
enter in the process only via the exchange interaction with s electrons on each ferromagnetic 
electrode. In the context of second quantization and neglecting the magnetization energy 
(Zeeman term), the unperturbed Hamiltonian is given by 

H = £ (1) 

k(T,a=(L,i?) 
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with L(R) referring to the left (right) ferromagnetic electrode, c££ (c£ CT ) are the creation 
(annihilation) fermionic operators for wave vector k and spin a, = 4^ — aA a is the 
Hartree-Fock energy and A a is the shift in energy due to exchange interaction in each side 
of the barrier. 

In writing the interaction part of the total Hamiltonian, which makes possible the transfer 
of electrons from one side of the insulating barrier to the other, we follow Ref. 8. Apart from 
the direct transfer which comes from elastic processes, we include transfer with magnetic 
excitations that originates from the s-d exchange between conduction electrons and localized 
spins at the interfaces. The excitations are described by a linearized Holstein-Primakoff 
transformation 16 , in the spirit of a one-magnon theory. We use the following Hamiltonian: 

Hint = E *ktf + + (2) 

kk'cr 

+7^ E *k q (<£[<£ T + c £t c y y^slb^ + v^j) + 

» s kk'q 

+7^ E ^k' q + ( + v^6?) + 

* s kk'q 

+^ E ^k'q (<£{<£, - cjt^ + fc. c .) (5 L + S R - (b%% + ftjtfcj)) , 

JVs kk'q 

where t^ k / is the direct transmission coefficient, is the inelastic transmission coefficient 
(depends on the exchange integral), S L (S R ) is the spin value at the left (right) side, N s is the 
total number of spins at the interface, and (b™) are the creation (annihilation) operators 
for magnons with wave vector q at each interface between the barrier and the electrodes. 
The wave vector q is quasi-two dimensional (the magnon wave function is localized at the 
interfaces, but with finite localization length). 

In general, the total current obtained with (2) has contributions from elastic processes, 
resulting in a direct tunneling which conserves spin, and from the inelastic ones, which 
involve emission and absorption of magnons with electronic spin flip. In the following we 
describe the direct term. 



III. DIRECT TUNNELING CURRENT: ANALYTICAL EXPRESSIONS 

Considering only the direct part of the tunneling process, which means elastic processes, 
without involving magnon excitations, the current is easily obtained by 2,3 ' 8,10 

ho = — J dE T d {E, V, d, $ ) W {c) (E,V)[f(E-eV)-f(E)} (3) 

where 

W {c) = Y.N*(E)N^E-eV) (4) 

and C denotes the configuration scheme, C = P for parallel and C = AP for anti-parallel, 
f(E) is the Fermi-Dirac distribution, N R and the density of states at the right and left 
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electrodes, respectively. T d (E, V, d, $ c 



ud 

\ L kk' 



is the tunneling coefficient, being a function 



of the energy E, the applied voltage V, the thickness of the barrier d and the barrier height 
$ - In fact, T d is a function of the overlap integral between the left and right wave functions 
inside the barrier region. 

The resistance is readily obtained by R = G^ 1 , where G = dl/dV is the differential 
conductance. In the low bias regime, we are interested in voltages smaller than the Fermi 
energy and only the states near the Fermi level will contribute to the transport, so we can 
expand the density of states in a Taylor series as follows: 



K(e) = E 



1 d n N*(E) 



n=0 



dE r > 



(E — Ei 



(5) 



Ep 



Now, let us calculate Wc for the P and AP configurations, using (5). In the P con- 
figuration the majority and minority bands in each electrode corresponds to the same spin 
orientation, and in AP configuration the majority band of one electrode is the minority on 
the other: 



(P) 



YY — 
^ ^ 



d*Nfi (E) dm^i (E - eV) 



dE* d{E- eV)' 



+ 



(6) 



+ 



<?N* (E) d>N^ (E - eV) 



dE i 



and 



W {A P) 



d(E-eV) 3 
d'N* (E) d?N^ (E - eV) 



(E - E F y (E-eV - E F ) 



+ 



dE* d(E- eVy 
d*N* (E) $N L M (E - eV) 



+ 



(7) 



d& d{E- eV) 3 



(E - E F y (E-eV- E F ) J 



E F 



Taking into account identical electrodes and the low bias regime, we can expand these 
expressions to first order with good accuracy. The s-band can be represented by a parabolic 
dispersion relation and density of states N a oc y/E — A a , where A <7 (<7 =t, I) gives the 
bottom of the spin band, with |Aj — A J = 2 A, as in Ref. 10. However, we consider here 
cases more general than the parabolic dispersion, with the band structure described through 
the following set of parameters: 




N 



M 



N m / F 
dN M /dE s 

dN m /dE / 

1 dN m \ 

N m dE j 



(8) 



with all quantities evaluated at the Fermi level, and m and M stand for minority and 
majority spin bands, respectively. We get the analytic expressions: 



W (P) = «) 2 { (l + r 2 ) + (3 (1 + r\) (2s -V)+p 2 (l + A 2 ) e(e - V)} , 



(9) 
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and 

W ( ap) = «) 2 {2r + (r + A) (2e - + /3 2 Ae( £ - V)} , (10) 

where e = E — E F and £ and V must be given in eV. 

There are several possibilities for including the tunneling transmission coefficient T d in 
the theory. One is the approach followed by Simmons 11 , where the barrier is parametrized 
by an effective height $o and an effective thickness d: 



T d (E, V, $ , d) = exp[--W2m($ - £*)] = exp 



■1.024 d\ & 



exp 



1 r/e 



(11) 



where all energies are measured from the Fermi level and given in eV, the barrier width given 
in Angstrom = 0.1 nm, and rj is some constant relating the energy e with its component 
e z perpendicular to the barrier. This latter parameter appears due to the fact that we are 
using a one dimensional formula to explain the behavior in the 3D case. 

Since the Fermi-Dirac functions are step-like at 0°K, we can easily obtain the conduc- 
tance for both configurations. 



G 



2vre 2 d 



(C) 



h dV \Jo 



V de T d (e,V,<S> ,d) W (c) (e,V) 



With some simplifications in the integration process (taking into account the behavior 
of the integrand in the range of integration, and making use of some geometric arguments), 
one obtains 

27re 2 
G(0 = — 



{A (c) T d (V,%,d) + 



(12) 



+ llv [ B (c)V 2 T d (V, $ , d) - C (c) V 3 T d (3V/5, $ , d) 



6dV 



being A(c),B(c) and C(c) constants related to the configuration scheme and the density 
of states. Following, the analytical expressions for the conductance in both parallel and 
anti-parallel configurations are presented, using (12) and considering the expansions (9) and 
(10): 



G 



2ire 2 



exp[-1.02W$ ] [Kf { (1 + r 2 ) exp[ 



h 



nVd 0(1 + r\) 



+ 



x 



X 



ndV 2 nVd nVd 
exp[^-^=J + l/(exp[ 



4v$o 2v$o 
(3 2 (1 + X 2 ) 



exp[ 



2^% 
3r]Vd n f jr2 



(13) 



10V$o 



V* + 



and 



G 



2ne 2 



(AP) 



h 



exp[-1.024dJ$ K] 2 { 2rexp[ 



r]V 3 d ' 



rjVd (3{r + A) 



x 



V dV 2 nVd f rjVd 

exp[ — j=\ + V exp[ 



4v$o 2v$o 

-/5 2 Aex P ; 



+ 



x 



lCK/C 



r]V 3 d 



(14) 
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The expressions above can be easily inverted to obtain the resistance, with the MR 
defined as 



Ai? Rap — Rp 



R 



Rai 



(15) 



Note that the above definition is limited to 100%, since Rap > Rp- In the limit V — > we 
have approximately 

G (P) = ^exp[-1.024d V / ^]ra 2 (l + r 2 )exp[^gL] 



and 



G t w, = ^exp[-1.024rf V / ^][<] 2 (2r)exp[^ 1 " 



h 



With the experimental value of AR/R at zero bias, we can easily obtain the ratio of the 
densities of states r at the Fermi level by: 



1 Ai? 




1 J? 





+ 



1 



1 



0-¥ 


) 2 




y=o/ 



- 1 



(16) 



which does not depend on the barrier parameters. In turn, the barrier height <3>o and thick- 
ness d determine the absolute value of the resistance. Typical values used in our examples 
are d = 1.0 nm and $ = 3.0 eV. Estimation of the resistance of such a junction yields 
resistance-area products of the order of RS ~ 3.3 x 10 4 [Q fim 2 ], where S is the junction 
area given in /im 2 . This value follows closely the resistance-area scaling obtained for differ- 
ent junctions in Ref. 17, with values of the MR ranging from 16% to 22 %. Representative 
experimental data of the tunneling resistance dependence on bias are given in Ref. 2, 3, 8, 
9. We compare our theoretical calculation with results presented in Ref. 8 at 4.2°K. There, 
the zero-bias MR is approximately of the order of 25%, which yields for the r parameter of 
(8) and (16) the value r = 2.21. In Fig. 1 we show our theoretical results for the resistance 
calculated with formulae (13) and (14). The band structure parameters were taken with the 
values A = 0.07 and f3 = 2.7, and the tunneling parameter as r] — 0.1. The small value of A 
depicts a situation where the majority spin band is saturated at the Fermi level, while the 
minority one has a large variation 10 . However, the slope of the resistance near zero bias only 
depends on the ratio of the densities of states, in the form 



Rap *Ro(l~ ±x) , 
Rp ~ Rq 



1 



1 



(17) 



where x 



rjd 



\V\ and Rq 



,1 + r 2 1 + r 2 
exp[1.024dv^o] 



T F]2 
1 ml 



is a scale factor related to the absolute 



V*o ' (¥) N> 

resistance. Note that we get a zero bias anomaly, but a good fit with the experiment is 
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only obtained for the parallel configuration, as in Ref. fO. One can adequate the theoretical 
model to a better fit with the data, using more terms in the Taylor expansion of Wm), 
or leaving the densities of states as free parameters 2 . However, we interpret the failure of 
fitting the data for the AP configuration as a hint that points to the contribution of an 
extra mechanism, which affects differently the P and AP resistances. The linear terms in 
(17) cancel out when one gets the MR, as shown in Fig. 2, along with the experimental 
data. We pursue our argument further in the next section, with the inclusion of magnon 
inelastic scattering processes in the calculation of the MR. 



IV. MAGNON- ASSISTED INELASTIC TUNNELING 

In this section we consider not only the elastic (spin conserving) processes but inelas- 
tic magnon- assisted contributions to the tunneling current. The latter are responsible for 
opening the spin-flip channels, substantially reducing the MR near zero bias. Magnons are 
spin-wave excitations 16 which interact with electrons, being emitted or absorbed, thus pro- 
ducing changes in their energy and allowing for spin-flip scattering. Electrons accelerated 
by the electric field relax their energy, producing those collective excitations at the magnetic 
electrode interfaces. At low temperature, only magnon emission processes give a significant 
contribution to the resistance. However we analyze in the following the general case, de- 
scribing each one of the eight processes associated with emission and absorption or magnons. 
There is one extra term related to the overlap between wave functions of the electrodes, not 
involving changes in the number of magnons. This term is proportional to the exchange 
transmission coefficient T J = \t{ k , q \ 2 , resulting in a very similar formula to the one found 
for the direct tunneling in the previous section: 



J de {S 2 R + S 2 L ) T J ( e , V, d, $ ) W {c) (e, V)[f(e -V)- f(e) 



Let us consider now the electron tunneling from the left to the right electrode with the 
emission of one magnon at the right side interface: 



2vre 

x (e -V + hcu)N r R (e) p7V)[l + f B E^)\f{e -V + hu)[l- /(e)] , 



Ui\ = — J duo J de2S R T J {e,V,d,<5> ) x 



being p ma9 (u) the density of magnons at the right side interface and Jbe the Bose-Einstein 
distribution: 

f 1 
lBE exp[^]-l 



An identical expression appears when considering the magnon emission at the left side 
interface yielding: 

I?c\ = ^jdujjde 2(^p7» + S LP T 9 ^)) T J (e, V, d, $ ) x (18) 



x 



Af(e -V + hu)N*(E) [1 + f BE (uj)]f(E -V + tuv)[l - f(e)\ . 
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When the tunneling occurs from right to left with one magnon emission, we have: 

Ifc] = 2 -^Jdujjde 2(S RP r r^) + S L pT 9 {^)) T J (e, V, d, $ ) x (19) 
x N{(e -V + hu)N?{e) [1 + f BE (u)]f(s)[l -f(e-V + hu)] . 

In turn, for magnon absorption we get: 

Itc) = 2 -^JdujJde 2{S RP T 9 ^) + S LP T 9 ^)) T J (e, V, d, $ ) x (20) 
x N L( £ _ V _ nw)N*(e) [f B E^)]f{e -V-hw)[l- f{e)} 

and 

1$) = 2 ^jdujjde 2(S R p r rH + S LP T 9 ^)) T J (e, V, d, $ ) x (21) 
x TVf ( e _ y _ ^)7Vf (e) [/ BB (w)]/(e)[l - /(e - F - Ml . 

The total current due to one magnon exchange is then: 

J m«9 _ 1 (C) + J (C) J (C) + J (C) _ J (C) 

Typical IET magnon spectra are shown by Ando and coworkers in Ref. 1. They display 
a strong peak around 12 — 20 mV and a rapid decrease for energies below the peak, due 
probably to a low energy cutoff, with a vanishing magnon density of states at very small 
energies. Introducing this low energy cutoff in the magnon spectrum, and taking the low 
temperature limit T — > 0°K, we get /be — ► for the Bose-Einstein distribution. This limit 
excludes the absorption terms in (22), leaving only the emission contributions to the total 
current: 

I mag = ^jdu J Q V ^dE{T J (e,V,d,%)(S R p ma9 (u) + S LP r 9 (io)) x 

xJVf(e - V - huj)Nf(e) (V - hco) - 

-T J (e,V,d^ )(S R pT 9 ^) + S L pr 9 ^)) x (23) 
xiVf (e)Nf(s — V — Tiuj) 6 (hu - V)} + ifa , 

being Q(x) the step function. 

One can use as the magnon dispersion relation a simple isotropic parabolic dependence, 
i.e., huj = E m (q/q m ) 2 , where E m is related to the Curie temperature by the mean field 
approximation E m = 3k R Tc/{S + 1), and q m is the radius of the first Brillouin zone 8 . In 
other words E m is the maximum magnon energy (high energy cutoff) 1 . Considering the 
above discussion, assuming identical ferromagnetic electrodes, and after some mathematical 
simplifications, one finally gets the conductance in the form: 

G c = G d c + G™ ag 

where G d c is given by (13) and (14), for P and AP alignment respectively, and G™ ag is shown 
below: 
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2vre 2 



(V) 2S 2 W (P) + A(V) SW {AP) 



(24) 



h 



and 




2vre 2 



(V) 2S 2 W (AP) + 2A(V)SW {P) 



(25) 



ft 



with 5 = (S R + S L ) /2 and 



A(V) 



V/E m iorV<E m , 
2 - E m /V for V > E, 



'in 



The functions W(p) and W^ap) in (24) and (25) have been evaluated near zero bias, using 
formulae (9) and (10) respectively, substituting e by the constant value 0.1 eV. We point out 
that magnon processes significantly contribute to the conductance only for voltages below 
100 mV, and so we can consider W(p) and W(ap) as almost constant under the integration 
sign. The exchange tunneling coefficient T J generally is one or two orders of magnitude 
smaller than the direct coefficient. We found excellent agreement between our theory and 
the experimental data using the same set of parameters of Fig. 1 for the tunneling barrier and 
the electronic structure, spin S — 3/2 and T d /T J = 37 for the ratio of the direct tunneling 
to the exchange tunneling coefficient. The magnon cutoff E m was taken to be 90 meV. 
The results are shown in Fig. 3 for the resistances and in Fig. 4 for the corresponding 
MR. Clearly, the AP configuration is more sensible to the magnon contribution, since the 
current for that configuration is weighted by the product N^N^j, which is much bigger 
than the factors N^N^ or N^N R which appear in the P current, with the indices m and 
M referring to minority and majority spin bands (compare Fig. 1 and Fig. 3). Obviously, 
minor differences between theory and experimental data come from the fact that we are 
using a very simplified model for the band structure and the magnon dispersion relation. 
We comment on these results in the next section. 



We have presented a consistent study of the voltage dependence of the 'giant' magne- 
toresistance in ferromagnetic tunneling junctions. Our approach includes: a) lowering of 
the effective barrier height with the applied voltage; b) different variations of the density of 
states for each spin band with voltage; and c) magnon assisted inelastic tunneling near zero 
bias. We found that taking into account all those effects is essential to fully explain exper- 
imental results at low temperature for the voltage range between and 500 mV. We have 
also clarify the role of the different parameters used in the theory: some of them (d, &o,v) 
determine the absolute value of the resistance at zero bias, which in turn is a scale factor 
in the theory; a different set, related to the band structure (r,/3,X), mainly monitors the 
global behavior with voltage and the value of the junction MR. To adjust our results with 
selected experimental data, we have taken f3, X > 0, but as shown in Ref. 10, this scenario 
is not unique and depends on the topology of the bands that contribute to the current; 
and finally, the behavior near zero bias (zero bias anomaly), with a rapid decrease of the 



V. CONCLUSIONS 
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resistance for the AP configuration up to 100 mV, is ascribed to magnon-assisted tunneling. 
Our estimation of T d /T J ~ 40 seems to be more realistic than previous estimations 8 . We 
note that the latter is the only adjustable parameter to fit the voltage dependence below 
100 mV (for both, P and AP configurations). Our calculation is in excellent agreement with 
the experimental data (see Fig. 3 and 4). 

Temperature effects are not discussed in this paper. As shown in Section IV, only magnon 
emission processes are included at low temperature (T — > 0). At finite temperature, we 
expect a decrease of the resistance near zero bias, due to one-magnon-absorption assisted 
tunneling 18 . The above should be superimposed to the thermal smearing in the Fermi-Dirac 
distribution of tunneling electrons 3 . 
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FIGURES 



Fig.l Resistance as a function of the voltage bias for the AP and P configurations: the 
experimental results (dotted line and symbols) are taken from Ref. 8 and the theoretical 
ones (solid lines) are calculated with formulae (13) and (14), using the following parameters: 
d = 1.0 nm, $o = 3.0 eV, = 1.0 in normalized units, r = 2.21, A = 0.07, j3 = 2.7, and 
7] = 0.1. The resistances are given in arbitrary units, normalized to the peak value at zero 
bias. 

Fig. 2 Magnetoresistance as a function of voltage. Parameters are kept the same as in 
Fig.l. 

Fig. 3 Resistance, in arbitrary units, as a function of the voltage bias for the AP and 
P configurations: the experimental results (dotted line and symbols) are taken from Ref. 
8 and the theoretical ones (solid lines) include magnon-assisted tunneling. Parameters are 
kept the same as in Fig. 1, with the addition of T d /T J = 37. 

Fig. 4 Magnetoresistance as a function of voltage. Parameters are kept the same as in 
Fig.3. 
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